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The Fulde-Ferrel-Larkin-Ovchinnikov (FFLO) state near the antiferromagnetic quantum crit- 
ical point (AFQCP) is investigated by analyzing the two dimensional Hubbard model on the 
basis of the fluctuation exchange (FLEX) approximation. The phase diagram against the mag- 
netic field and temperature is compared with that obtained in the BCS theory. We discuss 
the influences of the antiferromagnetic spin fiuctuation through the quasiparticle scattering, 
retardation effect, parity mixing and internal magnetic field. It is shown that the FFLO state is 
stable in the vicinity of AFQCP even though the quasiparticle scattering due to the spin fiuc- 
tuation is destructive to the FFLO state. The large positive slope dffppLo/dT and the convex 
curvature (d^i/ppLo/dT'^ > 0) are obtained, where //fplo is the critical magnetic field for the 
second order phase transition from the uniform BCS state to the FFLO state. These results 
are consistent with the experimental results in CeCoIns. The possible magnetic transition in 
the FFLO state is examined. 
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The BCS theory assumes the Cooper pairs which have 
zero total momentum.^ The Cooper pairing with finite 
total momentum was predicted in 1960's by Fulde and 
FerreP and Larkin and Ovchinnikov."^ The FFLO state 
can be stabilized in the presence of the high magnetic 
field or the molecular field. In addition to the U(l)-gauge 
symmetry the spatial symmetry is spontaneously bro- 
ken in this novel quantum condensed state. Although 
the FFLO state had not been experimentally observed 
in about 40 years after the theoretical prediction, sev- 
eral evidences for the FFLO state have been obtained in 
the recent experiments.* The FFLO superconductivity 
or superfluidity is attracting growing interests not only 
in the field of superconductivity* but also in the study 
of cold atom gases^ and high density quark matters.^ 

Vast studies of the FFLO state were triggered by 
the discovery of novel superconducting (SC) phases 
in the strongly correlated electron systems.* For ex- 
ample, a new SC phase has been identified in the 
heavy fermion superconductor CeCoIns at low temper- 
atures and high magnetic fields. ^""'^^ Experimental data 
for the FFLO state have been obtained in several or- 
ganic materials, such as A-(BETS)2FeCl4,i^ k-(BEDT- 
TTF)2Cu(NCS)2^^'^^ and (TMTSF)2C104.^°'2^ Lr these 
materials the paramagnetic effect which is essential for 
the FFLO superconductivity is important because the or- 
bital de-pairing effect^^ is suppressed by the heavy quasi- 
particle mass and/or the layered crystal structure. 

Another common characteristics in these compounds 
is the presence of the AFQCP in the vicinity of the SC 
phase. Although the critical spin fluctuation significantly 
affects the electronic state near the AFQCP, the influence 
of AFQCP on the FFLO state has not been investigated 
from the theoretical point of view. The main purpose 
of this paper is to examine the stability of the FFLO 
state near the AFQCP. Although the FFLO state has 
been investigated by Yokoyama et al.^^ on the basis of 
the random phase approximation (RPA) , the retardation 
effect, quasiparticle scattering and internal field, which 



play an essential role in our study, are neglected in the 
RPA. 

We analyze the two-dimensional Hubbard model in or- 
der to investigate the FFLO superconductivity near the 
AFQCP. The model is described as follows, 

H = XI ^(^)4^^fe<T + U Xl"'T"ii "hj^'^rii^- (1) 

We consider the square lattice and choose the following 
tight-binding dispersion, 

e(k) ~ — 2i(cos + cos ky) + At' cos fcx cos ky — /i, (2) 

in which the d-wave superconductivity is stabilized by 
the antiferromagnetic (AF) spin fluctuation. ^*"^^ In the 
following, we choose the unit of energy as 2t = 1 and fix 
t' /t = 0.25. The chemical potential is chosen so that the 
number density of electrons is ti = 0.9. The last term in 
eq. (1) describes the Zeeman splitting where h = ^q^bH 
with H being the magnetic field. As for the influence of 
the magnetic field, we here focus on the Zeeman split- 
ting which is essential for the FFLO superconductiv- 
ity.^'^ We assume a sufficiently large Maki parameter 
a = ^/2H°^'/H^2 and ignore the orbital effect^^' 27-33 f^^. 
simplicity. 

The model eq. (1) is analyzed on the basis of the FLEX 
approximation^* which is based on the Luttinger-Ward 
formalism. The critical fiuctuation near the AFQCP is 
taken into account in this approach. ^^'^^ xhe thermody- 
namic properties, such as the free energy, condensation 
energy and optical integral have been investigated at zero 
magnetic field. 

It is straightforward to extend the Luttinger-Ward for- 
malism to the FFLO state. We obtain the thermody- 
namic potential as follows, 

Vl{T,^i,h)^^Q{T,^i,h) + VLF + Q.B, (3) 

VIf = -X[log{detG(fc)-VdetG(°)(fc)-i} 

k 
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+ Y,G„{k)J:i{k) ~ F{k)A\k) - F\k)A{k)], (4) 

nB = ^G,,F,F%u (5) 

where iloiT, fi, h) = -T ^^^^ log[l+exp{- (3{e{k)-ah)}] 
is the thermodynamic potential at [/ = 0. We describe 
the normal and anomalous Green functions as, 



G(fc) 



_ f GT(fc+) 

F^k) 

_ ( G'f(fc+)-i 
At(fc) 



m 

-Gii-k-) 

A(fc) 

-G"(-fc_)-i 



T/NJ2n,, g- where uj„ = (2n + 1)t:T, n„ = 2m:T, T is 
the temperature and N is the number of sites. The unit 
h = c = kB = I is used throughout this paper. 

We carry out the numerical calculation on the N — 
64 X 64 sites and takes the Matsubara frequency A^f = 
2048 for T > 0.004. We keep the cutoff of Matsubara 
frequency as {N{ + l)7rT > 12.8 for T < 0.004, which is 
larger than the band width W = 4. 

(6) In the following results, the dominant component of 
SC order parameter A(fc) has the dx^-y^-wave symme- 
try. The order parameter has an admixed odd frequency 

(7) 2 -wave component in the finite magnetic field ow- 



where G^(fc)-i = gS°^(A:)-i - with Gi°^(fc) = 

(iwn — e{k) + ah)^^. We have described k± = k ± q^jl, 
and gp is the total momentum of Cooper pairs in the 
FFLO state. For simplicity we assume the Fulde-Ferrel 
(FF) state where the SC order parameter has the spa- 
tial modulation A(r) = exp(±igFr). Since we neglect 
the Larkin-Ovchinnikov (LO) state, the stability of the 
FFLO state is underestimated in the following results. 
Note that the FF state has the same critical magnetic 
field Fl^^ as the LO state when the second order phase 
transition is assumed. We will investigate the LO state 
in the vicinity of the AFQCP in another publication.^^ 

The normal and anomalous self-energies in eq. (7) are 
obtained from the generating function ^\Ga,F,F'^\ as, 

(5$ 



A(fc) 



8Ga(ky 



AT(fc) = - 



(5$ 



(9) 



(10) 



8F{ky 

by which the stationary conditions are satisfied as,"^^'^^ 

m m 



0. 



11 



SY.^{k) SA{k) (5At(fc) 

In the FLEX approximation the generating function is 
given as, 

$[G.,F,i^t] = Un^n^ +J2{iog[l ~ Uxliq)] 

q 

+ lu'[xUqr + X",iq)xliq)] + U[xl{q) + xUq)]}, (12) 
where the irreducible susceptibilities are described as, 
Xliq) - - Et^tCfc + q)Glik) + Fik + q)F\-k)], 



k 

X^(9) = -E^^(^ + 9)n^)- 



(13) 
(14) 

(15) 



ing to the violation of time reversal symmetry. The 
even frequency p-wave component is also admixed with 
the (ix2-y2-wave component in the FFLO state due to 
the breakdown of inversion symmetry like in the non- 
centrosymmetric superconductors.'^*''^® The extended s- 
wave component is also induced in the FFLO state unless 
qp II [110]. We will show that the p-wave order parameter 
plays an quantitatively important role for the stability of 
the FFLO state, while the other admixed components are 
negligible. In the following we assume the FFLO mod- 
ulation vector qp along the [110]-axis for the numerical 
accuracy, although a slightly larger condensation energy 
is obtained for qp \\ [100]. We have confirmed that quali- 
tatively the same results are obtained for qp \\ [100]. The 
direction of the modulation vector qp is actually affected 
by the orbital effect'^"' '^^ which is neglected in this paper. 
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We have used the notations J2k = S t ^^^'^ So ^ 
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Fig. 1. (Color onUne) Paramagnetic critical magnetic field 
for U/t = 2.8, 3.1 and 3.4 from the bottom to the top. The solid 
and dashed lines show the for the FFLO and BCS state, re- 
spectively. The tricritical point {T^pj^Q, ^^ppLo) shown by the 
circles (O)- The stars show the tricritical point {Tppj^Q, ^^^ppLo) 
obtained by neglecting the qp dependence of the order parame- 
ter A{fc). The Neel temperature Tn is shown by the dash-dotted 
line for U/t = 3.1 and 3.4. 



First we investigate the tricritical point 
(TppLQ, TJpplq) where the transition line from the 
FFLO state to the uniform BCS state merges to the 
critical magnetic field iJ^(T). Figure 1 shows the results 
of i?P(r) and (Tp*pLo,ifFFLo) ^r U/t = 2.8, 3.1 and 
3.4. We show the transition temperature Tn from the 
normal state to the AF state in the same figure.'"' 
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The Neel temperature Tn is determined by using the 
criterion 1 — Ux±{q) = 0.005 at q = (tt, tt) and fin ~ 0, 
where the AF order occurs with the magnetic moment 
perpendicular to the magnetic field. We find that the 
paramagnetic limit of the critical magnetic field 11^2 is 
significantly enhanced by the internal field arising from 
the AF spin fluctuation. This is consistent with the 
experimental results in CeCoIns.^^ 

We see that the FFLO state is stable in the vicin- 
ity of the AFQCP. The temperature at the tricritical 
point TppLQ in Fig. 1 (see also Fig. 2) is larger than 
that obtained in the weak coupling BCS theory Tpp^Q = 
0.56Tco^'''i where Tco is the critical temperature at zero 
magnetic field. This result seems to be surprising because 
the FFLO state is destabilized by the finite quasiparticle 
lifetime t* {k)^'~''^'^''^^ which is significantly decreased by 
the critical spin fluctuation. We here obtain the scatter- 
ing rate r,(fc) = 1/T*{k) = -z^{k)lm^^{k,0) - 2.3Teo 
at k ^ (7r,0) and a =T for U/t = 3.1 and iT,H) = 
(T^FLO^HhLo)- Here, z^^k) = {l-dRe^^{k,uj)/duj)-^ 
is the mass renormalization factor. This value of r(fc) is 
larger than the critical value F ^ 1.4Tco where the FFLO 
state disappears owing to the impurity scattering. 

The influence of the short lifetime is overcame by the 
other strong coupling effects, such as the retardation ef- 
fect and the parity mixing. We have confirmed that the 
frequency dependence of the SC order parameter A(fc) 
significantly increases the ratio Tpp^Q/Tco. The parity 
mixing between the d^a-yS-wave and p-wave order pa- 
rameters furthermore enhances the FFLO superconduc- 
tivity as shown in Fig. 1. The stars in Fig. 1 show the 
tricritical point (T'fflO' -^fflo) where the momentum 
and frequency dependences of the order parameter A(fc) 
is fixed to be those in the uniform BCS state. We see 
that the FFLO state is stabilized by the parity mixing, 
particularly near the AF critical point. Thus, the FFLO 
superconductivity is stabilized by the strong coupling ef- 
fects near the AFQCP. This is consistent with the fact 
that the experimental indications for the FFLO super- 
conductivity have been obtained in the strongly corre- 
lated electron systems near the AFQCP.**' '''^'^''"^^ 
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Fig. 2. (Color online) [/-dependence of the ratio TppLQ/TcO 
which represents the stability of the FFLO state. We show the 
Tppj^q/TcO for a comparison. 



Next, we discuss the J7-dependence of the FFLO in- 
stability. While the closeness to the AFQCP can be ex- 



perimentally controlled by the pressure or carrier doping, 
that is represented by the parameter U in our calculation. 
Figure 2 shows the ratio Tfplq/TcO which decreases as 
increasing U and approaching to the AF instability. This 
is mainly because of the quasiparticle scattering which 
is more sensitive to the closeness to the AFQCP than 
the other strong coupling effects. Thus, the FFLO state 
becomes unstable little by little as approaching to the 
AFQCP. This result is consistent with the pressure de- 
pendence of CeCoIns,*^ where TppLoZ-^cO gradually in- 
creases as increasing the pressure. 
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Fig. 3. (Color online) H-T phase diagram for U/t = 3.1 at high 
magnetic fields and low temperatures. Triangles show the second 
order phase transition from the uniform BCS state to the FFLO 
state. Diamonds show the second order phase transition to the 
AF state. We described the schematic thin dashed line below 
which the FFLO superconductivity coexists with the AF order. 
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Fig. 4. (Color online) H-T phase diagram of the BCS model 
eq. (16). The solid line shows the -ff^. The FFLO state is stabi- 
lized in the narrow region between the solid and dashed lines. 



Finally, we investigate the phase transition from the 
uniform BCS state to the FFLO state. Figure 3 shows 
the H-T phase diagram for U/t = 3.1. It is shown that 
the critical magnetic field i?FFLo(2^) from the uniform 
BCS state to the FFLO state rapidly decreases as de- 
creasing the temperature. The positive value of the slope 
di?FFLo/dT is much larger than that in the weak cou- 
phng BCS theory."*^ This is because the SC gap leads to 
the spin gap in the magnetic excitation and decreases the 
quasiparticle scattering rate. FFLO superconductivity is 
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enhanced below Tc since the Ufetime of quasiparticles 
significantly increases. We obtain the convex curvature 
(d^iypFLo/dT^ > 0) which is consistent with the exper- 
iments in CeCoIns.^'^^ The slope diJpFLo/dT decreases 
as decreasing the temperature because the quasiparti- 
cle scattering rate almost vanishes at low temperatures 
in the SC state. We note that the concave; curvature 
d^i?FFLo/dr^ < is obtained in Ref. 31 where the or- 
bital effect is taken into account but the influences of the 
spin fluctuation are neglected.'*^ Thus, the characteristics 
of the phase transition from the BCS state to the FFLO 
state in CeCoInj are reproduced by taking into account 
the critical spin fluctuation, while they seem to be incom- 
patible with the weak coupling BCS theories.'*' ^^'^^'^^ 

To illuminate the influence of the spin fluctuation more 
clearly we show the H-T phase diagram in the BCS the- 
ory. We obtain Fig. 4 by solving the following model, 



ariio 



fc.cr 



k,k' 



(16) 



on the basis of the BCS mean field theory. Here, B{k) = 
4'd{k)cj:^ ^c_g ^ and (f)d{k) = cos fcx — cos fcy. This model 
has been used in the weak coupling theories for the d- 
wave FFLO state. ^^'^^ We adopt the dispersion relation 
e{k) as eq. (2) and choose the attractive interaction Qs < 
so that Tco = 0.04. Wo assume the FF state and qp || 
[110] to compare Fig. 4 with Fig. 3 on an equal footing. 
We see that the FFLO state is stable in the narrow region 
just below Tc in Fig. 4. Thus, the FFLO state is much 
more stable in the repulsive Hubbard model (Fig. 3) than 
the BCS model (Fig. 4). This comparison illuminates the 
enhancement of the FFLO state near the AFQCP. 

The phase diagram in Fig. 3 seems to be consistent 
with the experimental results in CeCoIns, but there re- 
mains a discrepancy between our results and experi- 
ments. Although the first order phase transition from 
the normal to SC state occurs in CeCoIns above T = 
TppLQ,^^' the SC transition is in the second order in our 
calculation for T > 0.0032. This is because the quasipar- 
ticle scattering^* as well as the internal magnetic field^^ 
arising from the AF spin fluctuation suppress the first 
order phase transition. Therefore, the first order phase 
transition in CeCoIn5 above T = Tpp^Q should be at- 
tributed to the other origin rather than the AFQCP. For 
example, the orbital effect neglected in our calculation in- 
duces the first order phase transition above Tpp^Q.^^'^^ 
Furthermore, we found that the unequal effective mass 
between two spin species in CeCoIns^^'^^ also induces 
the first order phase transition, as will be shown in an- 
other publication.'^'' The heavy fermion physics, such as 
the strong electron correlation and the hybridyzation of 
localized and conduction electrons, may also affect the 
order of phase transition. 

At the last of this paper we comment on the magnetic 
instability in the FFLO state. Figure. 3 shows the phase 
transition from the FFLO state to the AF state (Dia- 
monds) in contrast to the result at the zero magnetic 
field where the AF spin correlation decreases as decreas- 



ing the temperature below Tc. Our result in Fig. 3 shows 
that the AF order can coexist with the FFLO state due 
to the large residual density of states in the FFLO state. 
We found that the AF order is furthermore enhanced 
in the LO state where the Andreev bound state at the 
spatial line node of SC order parameter induces the mag- 
netic moment.^'' These results may be consistent with the 
NMR measurement in CeCoIns,"'^ although the discrep- 
ancy seems to remain in the experimental data.**' *^'^^ 

In conclusion, we found that the FFLO superconduc- 
tivity is stable near the AFQCP even though the quasi- 
particle scattering rate is significantly increased by the 
critical spin fluctuation. The quasiparticle scattering is 
destructive to the FFLO state, but almost canceled by 
the strong coupling effect. The FFLO state is slightly 
suppressed as approaching to the AFQCP mainly be- 
cause the quasiparticle scattering is increased. The phase 
diagram obtained in our calculation is consistent with the 
experimental results in CeCoIns, however the first or- 
der phase transition is not reproduced. We pointed out 
some possible mechanisms for the first order phase tran- 
sition in CeCoIn5. The possible magnetic instability in 
the FFLO state was discussed. 

The authors are grateful to R. Ikeda, K. Izawa, Y. 
Matsuda for fruitful discussions. Nmncrical computation 
in this work was carried out at the Yukawa Institute 
Computer Facility. 
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